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We develop a model for the nonlinear oscillations of spherical drops composed of aqueous foam. 
Beginning with a simple mixture law, and utilizing a mass-conserving bubblc-in-ccll scheme, 
we obtain a Ray leigh-Plesset- like equation for the dynamics of bubbles in a foam mixture. 
The dispersion relation for sound waves in a bubbly liquid is then coupled with a normal 
modes expansion to derive expressions for the frequencies of eigenmodal oscillations. These 
eigenmodal (breathing plus higher-order shape modes) frequencies are elicited as a function of 
the void fraction of the foam. A Mathieu-like equation is obtained for the dynamics of the 
higher-order shape modes and their parametric coupling to the breathing mode. The proposed 
model is used to explain recently obtained experimental data. 

PACS numbers: 82.70.Rx,83.70.Hq 


I. INTRODUCTION 

Foams and froths are ubiquitous in nature and indus- 
try. They are the signature of vigorous, gas-entraining 
mixing processes in liquids. A minimalist conception of 
a foam would consist of a gas confined as bubbles within 
a liquid host. The aqueous foams considered here are 
composed of surfactant- bearing water and air bubbles. 
A comprehensive review of foam theories and applica- 
tions can be found in the textbooks by Edwards et al [1], 
Exerowa and Kruglyakov [2], and the article by Kraynik 
[3]. 

Theoretical treatments of the unique rheology of foams 
go back at least to Mallock [4], who was motivated to 
explain the common observation that “A tumbler con- 
taining a frothy liquid gives a dull sound when struck” . 
Mallock showed that the sound speed for intermediate 
void fractions was actually lower than its value for either 
the wet limit, c wa t er = 1500 m/s T or the dry limit, c a ir 
— 340 m/s. This result has been borne out by a century 
of subsequent work on bubbly liquids [5]- [9]. and, as 
seen later, leads to key insights into the free vibrations 
of foam drops. 
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The present work is motivated by the continuing need 
to measure, understand, model and eventually predict 
foam mechanics and rheology for wet or dry foams. To 
this end, two of the present authors recently described 
a non-contact technique in which small samples of foam 
(”foam drops”, see Fig. 1) were acoustically levitated 
and excited into resonance by modulating the levitation 
field [10]. Our technique utilizes acoustic levitation to 
provide both non-contact positioning and static and os- 
cillatory excitation of foam drops. By measuring the 
quadrupole eigenfrequency of a 3.8 mm radius foam drop 
to be 63 Hz, we inferred a shear modulus of roughly 73- 
78 Pa for a relatively dry foam. This value compared 
favorably with experimentally determined moduli utiliz- 
ing more traditional contact-based techniques ( [1 1]— [16] ) . 
The same technique has proven successful for determina- 
tion of the surface rheological properties of single-phase 
liquid drops[17]— [19] . 

The primary advantages of this acoustic levitation 
technique are its elimination of the requirement for sam- 
ple contact containment, its ability to test foams of arbi- 
trary gas volume fraction, and its ability to excite both 
shear and dilatational motion. The technique relies on 
a suitable physical model for the dynamic response of 
spheroidal foam drops. The heart of this model is the 
theoretical description of the foam material. Our earlier 
work [10] modelled the foam as an effective solid clas- 



tic medium. While successful in describing the small- 
amplitude oscillations of a dry foam with a fixed lattice 
of bubbles, such a material description has certain dis- 
advantages. First, a foam is only solid- like for high gas 
volume fractions and small amplitude motion. Second, 
such an effective elastic medium theory only implicitly 
incorporates the effect of varying gas volume fraction via 
the effective density of the medium, and it cannot capture 
the physics of a bubbly mixture. Finally, the gas volume 
fraction is a dynamic quantity, and during dilatational 
motion of significant amplitude it cannot be treated as 
a material constant. Wc thus need a model capable of 
describing wet foams. 

The subject of the present paper is a theoretical in- 
vestigation of the dynamics of wet foam drops. By us- 
ing a bubbly-fluid approach we obtain a model for the 
time-dependent response of a foam drop which retains 
the nonlinearity of the bubbles response. By lineariz- 
ing the equations we elucidate the resonance frequencies 
for the breathing and higher-order eigenmodes of a drop. 
We treat the coupling of shape and breathing modes, and 
investigate numerically the nonlinear equations obtained 
for the time evolution and interaction of the modal oscil- 
lations. A key feature is the inclusion of a time-varying 
void fraction in the equations. 

We begin in the next section with an argument for the 
relevant physics that must be included in the model. 

II. FOAM MECHANICS, RHEOLOGY, AND 
DROP DYNAMICS 

One of the most important characteristic parameters 
of a foam is its gas volume fraction ct g , or more commonly 
the 'void fraction/' A foam's thermodynamic, mechan- 
ical, acoustical and rheological properties are sensitive 
functions of the void fraction. Three regimes of foam 
morphology are typically identified. A “wet foam’ 1 (ap- 
proximately 0 < a g < 0.5) is essentially a bubbly liquid. 
The individual bubbles are free to move about within the 
liquid. Wet foams cannot support shearing motion, ex- 
cept at the surface of the individual bubbles. A “transi- 
tional” or ‘‘critical foam” (approximately 0.5 < a g < 0.7) 
is comprised of bubbles whose dynamics are strongly 
interacting, and whose surfaces may be in mechanical 
contact with each other. This regime may be usefully 
thought of as a phase transition between a liquid and 
solid-like state. A critical void fraction marks the point 
at which a foam begins to possess solid-like properties, 
such as shear wave propagation and yield stress. The crit- 
ical void fraction for three-dimensional foams is approxi- 
mately 0.67, which geometrically corresponds to random 
close packing of bubbles. Finally, a “dry foam 1 is the 
commonly encountered state in which the bubbles, at 
least for low to moderate straining rates, have a fixed 
position in a lattice. Such foams behave as viscoelastic 
solids for sufficiently small straining rates. However, a 


dry foam may (low as a liquid when strained beyond a 
critical point. 

Theoretical investigations of rheological properties be- 
gin with Derjaguin. [20] who derived an expression for the 
shear modulus of an idealized drv foam and showed that 
it was linearly proportional to the foam capillary pres- 
sure. This result implies that a foam’s modulus should 
scale as a tJ 1 *. Subsequent theoretical and numerical 
work has concentrated primarily on two-dimensional ge- 
ometric models limited to dry foams. The reader is re- 
ferred to references [21] [26] for examples of the ground- 
breaking work in this area and to reference [27] for a 
comprehensive review. Several theoretical models have 
addressed the unique rheological dependence on void 
fraction of foams. Bolton and Weaire [28] introduced 
a two-dimensional model that predicts the vanishing of 
the shear modulus at a critical void fraction. This model 
is strictly valid in the dry limit. In contrast, a static but 
bubble-based two-dimensional molecular dynamics sim- 
ulation [29] captures the transition features using wet- 
limit assumptions. The model assumes spherical bubbles 
that resist deformation because of their Laplace pres- 
sure. Among other things, the model predicts the vanish- 
ing of the shear modulus, but the scaling behavior near 
the transition is different than that, found for dry two- 
dimensional models. 

It is interesting to qualitatively consider the dynamics 
of foam drops in the limiting cases of wet, critical and 
dry. Wc consider a foam drop surrounded by a gas to 
simplify the situation. To engage in even a brief discus- 
sion. we must draw a distinction between the breathing or 
monopolo mode and the higher-order shape or multipole 
modes, since such motions are qualitatively different as 
well. First we consider breathing mode oscillations. The 
restoring force for perturbations from the drop's equilib- 
rium volume is provided bv the internal pressure of the 
individual bubbles within the drop, which will expand 
and contract when the drop volume is externally forced. 
Surface tension plays a small role, since the Laplace pres- 
sure 2cr / R q for a drop of radius Rq is much smaller than 
the ambient or atmospheric pressure. The mass is that of 
the liquid between the bubbles. Dissipation is provided 
by the bulk fluid motion, the surface fluid motions at the 
drop surface and at the individual bubble surfaces, and 
also by heat transfer and acoustic radiation of the indi- 
vidual oscillating bubbles [30]. For critical and dry foam 
drops, a primary difference is that surface tension be- 
comes more important as a restoring force because of the 
many thin film fluid connections which form inside the 
drop. The mass continues to decrease as the void frac- 
tion increases. It is difficult to make any general state- 
ment about the effect on dissipation due to increasing 
void fraction, except to say that dissipative effects are 
growing relative to inertial effects. 

For shape oscillations of wet foam drops, the restor- 
ing force for perturbations from the drop’s equilibrium 
shape is surface tension acting at: the drop interface. 
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Since surface- active agents are present, a local Marangoni 
restoring force due to gradients in surface tension also 
contributes. The effective mass is once again the mass of 
the liquid component. The dissipation is more strongly 
affected by the surface terms at the drop interface than 
for the monopole case, and the thermal and acoustic bub- 
ble dissipation terms are negligible. For shape oscilla- 
tions of critical and dry foam drops, the internal thin film 
fluid connections add stiffness to the drop, plus allowing 
the possibility of torsional multipolc modes. As for the 
monopolc case, the mass is decreasing, and dissipation is 
again ambiguous as the void fraction increases. 

The discussion above has implicitly assumed that the 
internal pressure of the mixture inside the drop is uni- 
form (except for the Laplace pressure contribution to 
the (again assumed uniform) interior bubble pressure). 
This will hold true until either the wavelength of incident 
sound is short compared to the drop radius R (which will 
never happen during standing wave acoustic levitation), 
or the velocity of the drop interface R approaches the 
speed of sound in the mixture, which is also unlikely. 
Uniform mixture pressure implies that all bubbles oscil- 
late essentially in phase unless there are wide disparities 
in the bubble size distribution, or nonlinear effects dom- 
inate. 

For both breathing and shape modes, foam drops near 
the critical void fraction will experience a sort of mode 
dispersion as bubbles begin to strongly interact with their 
nearest neighbors. Energy initially concentrated in a 
single global drop mode will be dispersed into motion 
of small collections of interacting bubbles. The result 
may well be an "apparent"’ increase in damping of the 
global observable drop modal oscillation. While this is 
highly conjectural, we believe we have indeed observed 
the apparent damping effect for forced quadrupole oscil- 
lations of near-critical foam drops. For the same liquid 
constituent, the quadrupole mode for both wet and dry 
foam drops was underdamped and thus a resonance was 
observable. However, for near critical void fractions the 
resonance was unobservable. A detailed investigation of 
this effect is a topic for future work. 

Thus we turn our attention to a more dynamic and 
bubble-based description of a foam, and in so doing we 
explicitly incorporate the fact that we utilize acoustics 
and acoustic levitation of bounded foam drops in our ex- 
periments. We wish to improve upon such existing mod- 
els by considering 3-dimensional cases, and by explicitly 
incorporating bubbly fluid dynamics and acoustic wave 
propagation in our model. We begin here by introducing 
a three-dimensional (spherically- or axisymmetric) model 
for the eigenmodal oscillations of spheroidal foam drops 
in the wet limit. 


III. WAVE EQUATION FOR WET FOAM 

In this section wc follow 8] to introduce the wave equa- 
tion for an aqueous foam in the wet limit.. For the en- 
visioned applications of this analysis, the liquid may be 
treated as incompressible. Viscous dissipation as it ap- 
pears in the normal stress balance for individual bubble 
oscillations is included. Thermal and acoustic dissipa- 
tion are neglected. Likewise, we do not consider processes 
such as bubble coalescence, breakup, or dissolut ion which 
affect the number of bubbles and/or their equilibrium 
size. 

To avoid confusion with the standard notation for ve- 
locity potentials (note that several authors use o to de- 
note the void fraction), let a./ and ck g be the fractional 
volume concentration, and p\ and p 9 the density of the 
liquid and the gas. respectively. Then the density p of 
the two-phase mixture is given by 

p = aipi - n 7 p 7 . (1) 

where ai + = 1. With p 9 <§C pi and the assumption 

that all bubbles have the same radius u, we obtain the 
following approximation for the mixture density: 

p « pi(\ - a,,), a,, = Ua V (2) 

where n is the number of bubbles per unit volume of the 
mixture. 

The dynamics of the bubbles in a foam may be treated 
with a sphcrically-symmctric bubblc-in-ccll scheme. Ac- 
cordingly, each single bubble in the bubbly liquid is 
placed in the center of a spherical cell with the radius 
A chosen in such a way as to cover the whole bubbly 
mixture by cells. 
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It should be noted that A.a,a 7 are variable in time and 
space, and a < A <C R typically. Here, R is the radius of 
the drop. 

Bernoulli’s integral for spherically-symmetric incom- 
pressible liquid motion in the ceil may be written in the 
following form: 

/ •• , 3 , Pi t 2-\ Pia*d 2 

P = Pa~ Pi\aa+ -a ) + — (a a) t (4; 

where dots denotes time derivatives, r is the radial co- 
ordinate with the origin at the bubble center, v is the 
shear kinematic viscosity of the liquid. p n is the liquid 
pressure at the bubble wall, which is relat'd to the gas 
bubble pressure p g by the formula 
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Pg — Pn 
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( 3 ) 


Here, a is the surface tension coefficient and v is the 
kinematic viscosity of the liquid. Often, it is treated as 
an “effective" viscosity which accounts for dissipation 8]. 

Taking into account that the liquid pressure at the 
cell boundary is equal to the pressure of the two-phase 
mixture p [5] one can obtain the following generalized 
Rayleigh-Plesset equation for radial bubble motion in the 
cell: 



( 6 ) 


The gas pressure in the bubble is assumed to be uni- 
form and variable in time and may be calculated via a 
polytropic approximation. 


Pg ~ 




(7) 


where k is the polytropic exponent (k — 7 ff for adia- 
batic and k = 1 for isothermal oscillations of the bubbles, 
where denotes the gas adiabatic exponent). 

Next we outline the approach to obtaining the mixture 
wave equation. The conservation equations for the num- 
ber of bubbles, the mass of the mixture, and momentum 
of the mixture are first linearized about the unperturbed 
state of the mixture. For this purpose we introduce small 
perturbations p\ p f . v* , n f . and a' of the equilibrium val- 
ues po, po = Pj( 1 - ago)- Fo = 0, no, and ao such that 
p — p 0 -(- p\ p = p 0 -j- p', u = v* , n = no + and 

a = a 0 + a! . Then we linearize equation (6) and omit 
any dissipative terms. Combining these Linearized bub- 
ble dynamics and mixture equations yields the (inviscid) 
equation of wave propagation in foam mixtures. Here 
and in the following we will write this equation in terms 
of a velocity potential p (v = Vp) keeping in mind that 
all variables p', p', and p are proportional to each other 
in the linear approximation. 
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Here cb is the speed of sound in a wet foam. *jsr is 
the frequency of free oscillations of single bubble in ail 
infinite Liquid. jJb is the frequency of free oscillations of 
bubbles in a foam. 


IV. NORMAL MODES OF A SPHERICAL 
FOAM DROP 

Let us consider a foam sample consisting of a spherical 
liquid drop of radius R with .V spherical gas bubbles dis- 
persed inside. To derive formulae for the cigcnfrequeneios 
of such a foam sample, let, the wave field inside the foam 
drop be described by 


p = exp(iu;£), (9) 

where uj is the (unknown) frequency of free oscillations. 
Substitution of (9) in the wave equation (8) results in 
the Helmholtz equation for the amplitude of the velocity 
potential 


V 2 <t> + A: 2 $ = 0. A J 



( 10 ) 


The oscillating surface of a foam drop can be expressed 
in the following way: 


oc 

r s = R o + ^2b n P n (cos 0) exp lujt, ( 1 1 ) 

n— 0 

where P n (cos0) are the Legendre polynomials, and the 
coefficients b n are the unknown amplitudes of eigenmodai 
shape oscillations ( b n Rq). 

The solution of Eqs (9). (10) 
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p = ^ B n j n (kr)P n (cos 0) exp iujt (12) 

n— 0 

describes the velocity potential inside the foam drop. 
Here j n are the spherical Bessel functions, and the B n 
are the unknown amplitudes for the acoustic field. 

The relationship between b n and B n can be found bv 
matching the normal component of the foam velocity at 
the drop wall with the normal displacement of the drop 
surface at that point, which in the linear case is 


dp 

dr 


= r*- 

r = /?o 


(13) 


From Eqs (12), (13) and (11) wc thus obtain a formula 
for the velocity potential in the foam drop 
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iub n jn(kr) 
k j'JkRo) 


P n (cos0)expL;f . 


(14) 


(6 ~ 10 -:{ ,,s ~ 10* '*). When h and s equal zero. Eq. (20} 
reduces to sin(r) — 0. for which the solution is r = m~. 
That means that Eq. (20) may be used to calculate a 
small correction to that solution. 


Bv imposing normal stress continuity we obtain the 
internal pressure 
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p = p 0 + — - pn 1 - a yO)-Q£ = 


(15) 


C = W7T+c. 5 <C 7T. m = 0, 1 (21} 

Substitution of Eq. (21) into Eq. (20) leads to the solu- 
tion 


~ — m~ 
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P n (cos0) exp iujt. 

( 16 ) 

That pressure must be balanced with sum of the external 
pressure and the Laplace pressure p „ 


2cr r-' {u — l)(n + 2) , 

Pa = — 4- <r > ^ 6 n P n (cos (9) exp iu;L 

Po 


n— 0 




(17) 

Equating terms with the same n we find the following 
equation: 


which gives the following generalized formula for the free 
monopole oscillations of a foam drop (compare to Eq. 
(24)) 
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(25) 

It is interesting to note that the effect of including surface 
tension is to lower the cigcnfrcquency. 

For Po large enough to ignore surface tension pressure 
(thus letting parameters b and s go to zero), we obtain 


2 _ &(n - 1 )(n + 2)kj , n (kRg) 
n P/(l-QV, 0 )Po Jn{kR oT 


( 18 ) 


Eq. (18) together with Eq. (10) constitute a full set of 
equations to calculate resonance frequencies of a foam 
drop. 

Let us consider two cases: monopolc (breathing) os- 
cillations (n = 0) and quadrupolc shape oscillations 
(n = 2). 


A Breathing mode oscillations 
When n = 0, Eq. (??) becomes: 


2 ffk j’o(fcflo) 
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(19) 
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Eq. (24) was evaluated for the parameters of a 
foam drop taken from the experimental observations [10*. 
Namely, P<} = 3.78 mm. 0.3 mm < ao < 0.5 mm. The 
normalized frequencies uj/u jsb versus foam void fraction 
a g o for m = 1 are shown in Fig. 1. The dashed line de- 
notes a gas volume fraction a g o = 0.77 of the foam drop 
used in [10]. One can sec that the monopolc frequency 
of the foam drop asymptotes to the single bubble result 
ujsb in the low void fraction limit, and to infinity for a 
void fraction of unity. In between these limiting values, 
the foam drop frequency is some fraction of ujsb- 

Typically, Po » ao s and Eq. (24) may be simplified as 
follows 


Using an explicit expression for the spherical Bessel 
function of zero order (jo(kRo) = sm(kRo)/kRo) this 
equation may be presented in dimensionless form as 


z 2 sin z 
1 -f bz 2 z 
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Thus, with the Eq. (8) and taking into account that 
2cr/ao <£ po one can get the formula for the frequency of 
the foam drop monopole oscillations, which looks similar 
to the well-known Minnaert formula [31] for single bubble 
monopolc oscillations 


For typical parameter values used in our experiments 
[10] the dimensionless variables b and s are very small 
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Here a* is an effective bubble radius. It is easy to es- 
timate that for the experimental data (Ro = -1.78 mm. 
a, v o = 0.77, 0.3 mm < a 0 < 0.5 mm) the effective bubble 
radius a* = 0.88 mm is larger than bubble radius and 
smaller than radius of the foam drop (a<) < a, <*Rq). 

These results, especially Eq. (24) may be usefully com- 
pared to previous results for the breathing mode of a 
compact bubble cloud in water derived by several au- 
thors [32: [35]. In those works, the motivation was to 
explain low frequency ambient noise in the ocean as due 
to collective oscillations of clouds of bubbles. The low 
void fraction limit is the same, but the high void frac- 
tion limit is not, since in the present case the high void 
fraction limit corresponds to the effective mass of the os- 
cillator approaching zero, whereas in the bubble cloud 
scenario, the high void fraction limit corresponds to a 
single bubble of radius R. 

It is remarkable that (approximately) the frequency 
does not depend on bubble radius, but depends only on 
the void fraction of the foam and the radius of the foam 
drop. It is clear that the apparent difference between 
frequencies for different bubble radii shown in Fig. I is 
only due to differences in the single bubble frequency ujsd 
by which we normalize. For instance, when a g o = 0.77 
(dash line) the frequencies of monopole oscillations for 
ao — 0.3 mm and no = 0.5 mm arc equal to 3708 Hz and 
3674 Hz respectively. This strengthens the case for our 
working assumption that the most important parameter 
is the void fraction - the details of the bubble size distri- 
bution should not affect the leading order results. 


B Shape mode oscillations 

When n = 2 Eq. (??) becomes: 
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This equation may be presented in dimensionless form as 
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where dimensionless variables z, b, s are defined in Eq. 

(20). 

To look for the low-frequency solution of this equation 
one should use the asymptotic expansion for the spherical 
Bessel function of second order when z — ► 0 [36], 
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Then Eq. (28) may be rewritten in the following form 


1 + bz 2 




There is one low-froquencv solution of this equation 
which is approximately equal lo 
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Here u)L 2 is the Lamb frequency of quadrupolo oscilla- 
tions for an inviscid and incompressible liquid drop [37] 
of density pi( 1 — a„o) surrounded by a vacuum. In our 
case u>L2 ~ 210.s“ l . Eq. (32) presents a correction to 
that frequency due to foam compressibility. For our ex- 
perimental parameters the correction term is negligibly 
small (~ 10“ l ). 


V. COUPLING OF BREATHING AND SHAPE 
OSCILLATIONS 

Here we consider the influence of breathing 
(‘monopole’ or ‘radial 1 ) oscillations on the quadrupole 
mode. This corresponds to an experimental situation 
in which the breathing mode is directly excited by an 
external field with pressure gradients large compared to 
the drop radius R. 

It is instructive to compare the period of quadrupole 
oscillations with the time required for an acoustic wave 
to cross the foam drop. At frequencies about 200 ,s “ 1 (see 
cj £,2 in the previous section) the period of oscillations is 
about 3 - 10“ 2 ,s. The speed of sound for our foam is 
cq ~ 30m/<s. Then for a foam drop with radius R. = 
3.78 mm, the acoustic transit time is about 10“ 4 s which 
is much less than the period of the quadrupole mode. 
Thus, one can assume the foam drop liquid pressure p ( i 
to be uniform in space and variable in time, as was done 
for individual bubbles in section 1 above. 

As in Section 3, each single bubble in the bubbly liquid 
is placed in the center of spherical cell with the radius A 
chosen in such a way as to cover the whole bubbly mix- 
ture by cells. The only difference is that in Section 3 all 
variables were functions of time and space coordinates. 
In the present case, all these variables depend only on 
time. Thus, Eqs (2), (3). (5), (7) are valid, but, for tem- 
poral variations only. 

Taking into account that the pressure of the two-phase 
mixture pd is uniform in space and time-dependent, one 
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(30) 


can get the following generalized Rayleigh- Plesset. equa- 
tion for radial bubble motion in a cell (compare to Eq. 
16)): 



Conservation of the (incompressible) liquid volume in 
the foam drop leads to the following simple relationship 
between bubble and drop radii. 


ft 3 - flg = N (a 3 - ag) • (35) 

Eqs. (3),(5)-(7),(refR-P-l), and (35) form a complete 
nonlinear system that describes radial motion due to 
foam drop oscillations. 

In case of weak (linear) oscillations the system may be 
reduced to a forced linear oscillator: 


The velocity distribul ion and potential inside the foam 
drop are: 

w = tt = + d: = ~ r ~ ^ nB f) (t)r Tl ~ l P„(ca*fl) 

or Ft 

n — 1 

(3!)) 

r’ = + 'P = -yjf + Un(t)r n P n (cos 9) (40) 

n= L 

Here B n (t) arc the unknown amplitudes of the foam mix- 
ture motion due to shape oscillatious. 

The relationship between b n (t) and B n (t) can be found 
by matching the normal component of the foam velocity 
at the drop wall with the normal displacement of the drop 
surface at that point, which in the linear approximation 
is 


R 4- JqR = — Fsincjt, (36) 


F - 

Pi (l - O!q0 3 ) a 0 

Here R = {R - Rq)/Rq is the relative foam drop radius 
fluctuation, ujq is the frequency of free oscillations of 
bubbles in foam (see Eq. (8)), Ap is an external driving 
(Pd = po + Apsinwt). The Laplace pressure on the drop 
wall is neglected. 

In order to investigate coupling between breathing 
and shape oscillations it is useful to consider the radial 
breathing motion as the base motion (below marked by 
an asterisk) and shape oscillations as a perturbation of 
this base motion (below marked by a tilde). In the homo- 
baric approximation discussed above, the velocity distri- 
bution tc*(r, t) for monopole oscillations (which satisfies 
the mass balance equation for the foam mixture) is given 
by 


^!r--r s = (-11) 

From Eqs (38), (39) and (41) we thus obtain the for- 
mulae for the velocity and velocity potential in the foam 
drop 


R r _ 
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In order to derive an approximate equation for drop 
shape oscillations let us consider the momentum equation 
for the foam mixture, 


p— r Vp-0. (44) 


w* = 


R 

~R' 


(37) 


The dynamic drop shape is expressed by the expansion 


where d/dt is the material derivative, and v is the veloc- 
ity of the mixture (liquid and bubbles move w r ith equal 
velocities). We now take a linearized form of the radial 
component of (44) as follows. 


oc 

r, = R(t) + J2 b n(t)P n (cos9), (38) 

n — 1 


w'hcrc P n ,(cos0) are the Legendre polynomials, R is the 
instantaneous foam drop radius and the amplitudes of 
eigenmodal shape oscillations b n (b n Rq) are unknown 
functions of time. Our goal is to derive a coupled set of 
differential equations for the variables R(t),b n (t). 


( i dw dw 

1 W + Wm lfr 




(45) 


Here p* = pi(l — a g ) is a function of time. Integration 
over the radial space coordinate r leads to the following- 
specific form of Bernoulli's integral 


P* 


li 



P = $(t). 


( 1 «) 


s 


The function ^(7) may he determined by evaluating the VI. NUMERICAL RESULTS 

left-hand-side of Eq. (46) at the limits r = r s and r = 0. 

We obtain approximately A Nonlinear breathing dynamics 


dt 


r — r ■ 


V ELE 

n dt 




P„( cosfl). 


TDD 


B Parametric instability 


c)\p 

dr 



P n (cos 0 ), 


W M \ 


r~r A 


R , 


Pir=r, * * f; (n ~ 1 £ n + 2) b n P n ( COS6). 


Then Eq. (??) leads to the following equation for the 
evolution of the amplitudes of the axisyminetric modes 


bn 4 - 


(n — l)n(n -h 2)rr 

pLr 3 



(47) 


which taking into account that p«R 3 = pi(i — a 3 o)i?o 
may be rewritten as follows 


bn 4 



(n — l)n(n H- 2)er 

(48) 


where LOL n is the Lamb frequency of shape oscillations of 
a foam drop with mode number n and effective density 
Pi(l — a 0 o). This result is analogous to the well-studied 
problem of parametric shape oscillations of single bubbles 
in an infinite fluid, see for example [38] and [39]. 

To obtian an analytical expression, if we consider only 
linear foam drop monopoie oscillations, then the solution 
of Eq. (36) is 


R = - 


Fsinu;£ 


- v* 


(49) 


It is beyond the scope' of this work to investigate dissi- 
pative effects. Nevertheless in the low void fraction limit 
we may incorporate the effects of weak liquid viscosity 
into the equation of motion for the shape oscillations de- 
scribed by Eq. (48). We give here without derivation the 
parametric equation of motion including weak damping 
of the shape oscillations: 



(n - 1 )(2n 4 1 

Pi ( I - o ) R{ 


where 3i n is the classical expression for the weak vis- 
cous damping of the shape oscillations of a pure liquid 
drop . We integrate Eq. (51) simultaneously with Eqs 
(34).(5)-(7),(35) to investigate the effect of a time-varying 
void fraction on the stability of the spherical foam drop. 

The regions of parameter space {K.e) where one ob- 
tains stable or unstable parametric shape oscillations are 
shown on Fig. 2, where stable zones are colored in black. 
One can see that the optimal driving frequency leading to 
parametric instability is two times larger than the eigen- 
frequcncy of the shape oscillations ( K = 1). the 2:1 ex- 
ternal resonance. 

In Fig. 3 we plot the frequency of quadrupolc oscilla- 
tions ijji 2 / 2 ir and the parametric driving resonance fre- 
quency jj — (jj r — 2u} l 2 versus foam void fraction 
For our experiments, the ^v 2 ) ^ 210 s“ l , the frequency 
of parametric resonance is jj r ^420 s “ 1 . It corresponds 
to a frequency f = ^65 Hz. which was observed 

in experiments [10]. 


VII. DISCUSSION 


and Eq. (48) may be presented in the form of Mathieu s 
equation. 


— r + ( K - 2e sin 2 r)6„ = 0, (50) 

dr * 


4 oj 2 r 

K = ~a' s ' 


2 F 


ut 

ji> — u; 2 ’ 2 


Solutions of Eq. (50) arc well studied (see [36]). There 
exist two types of solutions in the inviscid limit: bounded 
(stable) and unbounded (unstable) in time. 


A new bubble- based mathematical model of foam drop 
dynamics is presented. According to this model different 
types of foam drop oscillations are analyzed. The derived 
formulae arc used to explain the experimental data pub- 
lished recently [10]. In that work, an experimental tech- 
nique was demonstrated for acoustically levitating aque- 
ous foam drops and acoustically exiting their spheroidal 
modes. Results were presented in which a foam drop with 
gas volume fraction a 9 o = 0.77 was levitated at 30 KHz 
and excited into a quadrupole resonance at 63 ±3 Hz. 

Here it is shown that monopole (radial) oscillations 
and quadrupole (shape) oscillations, for foam parameters 
taken from experimental observations, are characterized 
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by two different scales of frequencies: about 3700 Hz for 
monopole oscillations, and about 33 Hz for quadrupole 
oscillations. 

Monopole frequency, being far from both, levitation 
frequency and quadrupole frequency, probably does not 
play any role in the phenomenon. The quadrupole os- 
cillation may be exited bv parametric resonance mech- 
anism. The highest frequency of the parametric exita- 
t ion is about 65 Hz. which corresponds with experimental 
data [101 very well. 
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